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We show that the traditional concept of the uniform electron gas (UEG) — a homogeneous system of finite 
density, consisting of an infinite number of electrons in an infinite volume — is inadequate to model the UEGs 
that arise in finite systems. We argue that, in general, a UEG is characterized by at least two parameters, viz. the 
usual one-electron density parameter p and a new two-electron parameter tj. We outline a systematic strategy to 
determine a new density functional E(p, tj) across the spectrum of possible p and tj values. 
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I. INTRODUCTION 

The year 2012 is notable for both the journal and one of us 
for, in the months ahead, both TCA and PMWG will achieve 
their half-centuries. It is inevitable and desirable that such 
occasions lead to retrospection, for it is often by looking back- 
wards that we can perceive most clearly the way ahead. Thus, 
as we ruminate on the things that we ought not to have done, 
we also dream of the things that we ought now to do. 

The final decade of the 20th century witnessed a major rev- 
olution in quantum chemistry, as the subject progressed from 
an esoteric instrument of an erudite cognoscenti to a common- 
place tool of the chemical proletariat. The fuel for this revo- 
lution was the advent of density functional theory (DFT) |T|| 
models and software that were sufficiently accurate and user- 
friendly to save the experimental chemist some time. These 
days, DFT so dominates the popular perception of molecular 
orbital calculations that many non-specialists now regard the 
two as synonymous. 

In principle, DFT is founded in the Hohenberg-Kohn theo- 
rem |2) but, in practice, much of its success can be traced to 
the similarity between the electron density in a molecule and 
the electron density in a hypothetical substance known as the 
uniform electron gas (UEG) or jellium ||3l 4T8l . The idea — the 
Local Density Approximation (LDA) — is attractively simple: 
if we know the properties of jellium, we can understand the 
electron cloud in a molecule by dividing it into tiny chunks of 
density and treating each as a piece of jellium. 

The good news is that the properties of jellium are known 
from near-exact Quantum Monte Carlo calculations IT9TI291 . 
Such calculations are possible because jellium is characterized 
by just a single parameter p, the electron density. 

The bad news is that jellium has an infinite number of elec- 
trons in an infinite volume and this unboundedness renders it, 
in some respects, a poor model for the electrons in molecules. 
Indeed, the simple LDA described above predicts bond ener- 
gies that are much too large and this led many chemists in the 
1970s to dismiss DFT as a quantitatively worthless theory. 

Most of the progress since those dark days has resulted from 
concocting ingenious corrections for jellium's deficiencies. 
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For example, significant improvements in accuracy can be 
achieved by using both the density p (r) and its gradient Vp (r), 
an approach called gradient-corrected DFT l30l . Even better 
results can be achieved by including a fraction of Hartree- 
Fock exchange (yielding hybrid methods OTl [321 ) or higher 
derivatives of p(r) (leading to meta-GGAs fl33"ll ). 

However, notwithstanding the impressive progress since the 
1970s, modern DFT approximations still exhibit fundamen- 
tal deficiencies in large systems l34l . conjugated molecules 
ll35l . charge-transfer excited states [36], dispersion-stabilized 
systems l37ll . systems with fractional spin or charge J38), 
isodesmic reactions [39| and elsewhere. Because DFT is in 
principle an exact theory, many of these problems can be traced 
ultimately to the use of jellium as a reference system and the 
ad hoc corrections that its use subsequently necessitates. It is 
not a good idea to build one's house on sand! 

In an attempt to avoid some of the weaknesses of jellium- 
based DFT, we have invented and explored an alternative 
paradigm called Intracule Functional Theory (IFT) 140-1421. In 
this approach, the one-electron density p (r) is abandoned in 
favour of two-electron variables (such as the interelectronic 
distance ry£) and we have discovered that the latter offer an 
efficient and accurate route to the calculation of molecular 
energies [43-48 1. Nonetheless, IFT is not perfect and has 
shortcomings that are complementary to those of DFT. As a 
result, one should seek to combine the best features of each, 
to obtain an approach superior to both. That is the goal of 
the present work and we will use atomic units throughout this 
article. 



II. ELECTRONS ON SPHERES 

In recent research, we were led to consider the behavior of 
electrons that are confined to the surface of a ball. This work 
yielded a number of unexpected discoveries ||49H57ll but the 
one of relevance here is that such systems provide a beautiful 
new family of uniform electron gases (see also [58 1). 

A. Spherium atoms 

The surface of a three-dimensional ball is called a 2-sphere 
(for it is two-dimensional) and its free-particle orbitals (Table 
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TABLE I. The lowest free-particle orbitals on a 2-sphere 
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TABLE II. Number of electrons in L-spherium and L-glomium atoms 
L 012345 67 



L-spherium 2 8 18 32 50 72 98 128 
L-glomium 2 10 28 60 110 182 280 408 



TABLE III. The lowest free-particle orbitals on a glome (i.e. a 3- 
sphere) 
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[Ijl are the spherical harmonics Y/ m (0, 0). It is known that 

2/ + 1 
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and doubly occupying all the orbitals with 0<1<L thus yields 
a uniform electron gas. We call this system L-spherium and 
will compare it to fwo-dimensional jellium [ 18 1. 
The number of electrons (Table |H} in L-spherium is 

» = 2(L+1) 2 (2) 

the volume of a 2-sphere is V = 4kR 2 and, therefore, 
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B. Glomium atoms 
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The surface of a four-dimensional ball is a 3 -sphere (or 
"glome" ED) and its free-particle orbitals (Table |llT|> are the 
hyperspherical harmonics Y[„ m (x,9,(j)). It is known 11601 that 
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and doubly occupying all the orbitals with 0<1<L thus yields 
a uniform electron gas. We call this system L-glomium and 
will compare it to f/iree-dimensional jellium lfT8l . 
The number of electrons (Table [IT| in L-glomium is 

« = (L+l)(L + 2)(2L + 3)/3 (5) 

the volume of a 3-sphere is V = 2k 2 R? and, therefore, 
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C. Exactly solvable systems 

One of the most exciting features of the two-electron atoms 
0-spherium and 0-glomium is that, for certain values of the 
radius R, their Schrodinger equations are exactly solvable (51]. 
The basic theory is as follows. 

The Hamiltonian for two electrons on a sphere is 
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where u is the interelectronic distance r\i = |i"i — r2 1- If we 
assume that the Hamiltonian possesses eigenfunctions that 
depend only on u, it becomes 
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where %> is the dimensionality of the sphere. Three years ago, 
we discovered that H has polynomial eigenfunctions, but only 
for particular values of R. (This is analogous to the discovery 
that hookium[? ] has closed-form wavefunctions, but only for 
particular harmonic force constants |62 63 1.) We showed that 
there exist |_(# + 1 ) /2J nth-degree polynomials of this type and 
that the associated energies and radii satisfy 



4R 2 n R n> m = n(n + 2@-2) 



(9) 



where the index m = l,...,[(n+ 1)/2J . 

For 0-spherium (i.e. S — 2), by introducing x — uj (2R) and 
using Eq. |9]l, we obtain the Sturm-Liouville equation 
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The eigenradii R can then be found by diagonalization in a 
polynomial basis which is orthogonal on [0, 1]. The shifted 
Legendre polynomials are ideal for this |64|. 
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For 0-glomium (i.e. 



d 
dx 



x(l-x 2 ) , s d>¥ 

w(x) — 

2 dx 



3), proceeding similarly yields 
- n ( n + 4 ^ X w(x)y=Rw{x)y (11) 



where the weight function w(x) = x\/\ —x 2 . The eigenradii 
are found by diagonalization in a basis which is orthogonal 
with respect to w{x) on [0,1]. 

An exact energy can be partitioned into its kinetic part 



£ T = (-1/4)(^|V2 + V2|vp) 
and its two-electron part 

£ee = (l/2)( V P|«- 1 | V P) 
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and the resulting reduced energies (i.e. the energy per electron) 
of the ground states of 0-spherium and 0-glomium, for the first 



two eigenradii, are shown in the left half of Table IV 



Many formulae have been proposed for Ex and E^ s but 
the most famous are those explicitly designed to be exact for 
^-jellium. In the case of exchange, one finds 
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where Dirac [5| determined the coefficient Xj, in 1930 and 
Glasser [68] found the general formula for X@ in 1983. 
The correlation functional is not known exactly but accu- 
rate Quantum Monte Carlo calculations on jellium in 2D [20j- 

[23][25] [26l[28]|29l and 3D GHElEl have been fitted 112611691 
to functions of the form 



„jell 
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By construction, Eq. (21 1 yields the correct energy when ap- 



plied to the uniform electron gas in jellium. But what happens 
when we apply it to a uniform electron gas on a sphere? 



III. SINGLE-DETERMINANT METHODS 



IV. THE NON-UNIQUENESS PROBLEM 



A. Hartree-Fock theory 



In the Hartree-Fock (HF) partition, the reduced energy[? ] 
of an H-electron system is 
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where the five contributions are the non-interacting-kinetic, ex- 
ternal, Hartree, exchange and correlation energies, respectively. 
The first four of these are defined by 
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where y/, (r) is an occupied orbital, p(r) is the electron density, 
and v(r) is the external potential. The correlation energy Eq is 



defined so that Eq. ( 14 1 is exact 



B. Kohn-Sham density functional theory 1 67] 



In the Kohn-Sham (KS) partition, the energy is 



£ K s = T s +E v +Ej+E x +E£ 



KS 
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where the last two terms, which are sometimes combined, 
are the Kohn-Sham exchange and correlation energies. The 
correlation energy E^ s is defined so that Eq. ( 19 1 is exact. 



The deeply disturbing aspect of jellium-based DFT models 
— and the launching-pad for the remainder of this paper — is 
the countercultural claim, that 

77ie uniform electron gas with density p is not unique. 

Though it may seem heretical to someone who has worked 
with jellium for many years, or to someone who suspects that 
the claim violates the Hohenberg-Kohn theorem, we claim 
that two ^-dimensional uniform electron gases with the same 
density parameter p may have different energies. To illustrate 
this, we now show that density functionals [5 26 69 1 which 
are exact for jellium are wrong for 0-spherium and 0-glomium. 



A. Illustrations from exactly solvable systems 

The energy contributions for 0-spherium and 0-glomium are 
easy to find. There is no external potential, so E\ = 0. The 
density p(r) is constant, so the Kohn-Sham orbital y/(r) = 
•v/p (r) is constant, and 7s = 0. The Hartree energy is the 
self-repulsion of a uniform spherical shell of charge of radius 
R and one finds 1 50 1 



T(2>/2 + \/2) 1 
R 



T(@-\/2) T(3>/2) 
The exchange energy is predicted [57| by Eq. (20 1 to be 
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and the correlation energy predicted by Eq. ( 2 1 1 is simply 
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(24) 



Applying these formulae to the exactly solvable states of 0- 



spherium and 0-glomium considered in Section II C yields 
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TABLE IV. Exact and Kohn-Sham reduced energies of the ground states of O-spherium and 0-glomium for various eigenradii R 



Exact Jellium-based Kohn-Sham DFT Error 





2R Ej E ee 


E 




Ej -E X -4 e " 


Eks 


Eks -E 


O-spherium 


V3 0.051982 0.448018 


1/2 





1.154701 0.490070 0.1028 


0.562 


0.062 




V2S 0.018594 0.124263 


1/7 





0.377964 0.160413 0.0593 


0.158 


0.015 


0-glomium 


vTO 0.014213 0.235787 


1/4 





0.536845 0.217762 0.0437 


0.275 


0.025 




V66 0.007772 0.083137 


1/11 





0.208967 0.084764 0.0270 0.097 


0.006 



the results in the right half of Table IV In all cases, the KS- 
DFT energies are too high by 10 - 20%, indicating that the 
correlation functional that is exact for the uniform electron gas 
in jellium grossly underestimates the correlation energy of the 
uniform electron gases in 0-spherium and 0-glomium. 



B. Limitations of the one-electron density parameter p 



The results in Table IV demonstrate conclusively that not all 
uniform electron gases with the density p are equivalent. The 
simplest example of this is the ground state of two electrons on 
a 2-sphere with R = \/3/2. The exact wavefunction, reduced 
energy and density of this system are 



¥ = 1 + u 
E = l/2 
p(r)=2/(37T) 



(25) 
(26) 
(27) 



but, when fed this uniform density, the exchange-correlation 
functional that is exact for two-dimensional jellium grossly 
overestimates the energy, yielding £ks = 0.562. 

This discovery has worrying chemical implications. Con- 
trary to the widespread belief that the LDA (e.g. the S-VWN 
functional) is accurate when applied to regions of a molecule 
where p (r) is almost uniform (such as near bond midpoints), 
our results reveal that it actually performs rather poorly. 

The discovery also has counterintuitive implications at a 
theoretical level. It implies that the years of effort that have 
been expended in calculating the properties of jellium do not 
provide us with a complete picture of homogeneous electron 
gases. On the contrary, although they inform us in detail about 
the infinite uniform electron gas, they tell us very little about 
the properties of finite electron gases. 

In a nutshell, the results in Table|IV|reveal that a UEG is not 
completely characterized by its one-electron density parameter 
p. Evidently, something else is required. . . 



C. Virtues of two-electron density parameters 

We know that it is possible for two uniform electron gases 
to have the same density p but different reduced energies E. 
But how can this be, given that the probability of finding an 
electron in a given volume is identical in the two systems? The 




FIG. 1. Specific Coulomb holes for 0-spherium (dotted) with R - 
\/3/2 and 2D jellium (solid). Both are uniform gases with p : 
2/(3w). 



key insight is that the probability of finding two electrons in 
that volume is different. 

This is illustrated in Figs [T] and |2j which compare the proba- 
bility distributions of the interelectronic distance [52 701I7TI 
in various two-dimensional uniform electron gases. These 
reveal that, although similar for u « (because of the Kato 
cusp condition 11721 ). the specific Coulomb holes (i.e. the holes 
per unit volume [40 1) in two gases with the same one-electron 
density p can be strikingly different. In each case, the jellium 
hole is both deeper and wider than the corresponding spherium 
hole, indicating that the jellium electrons exclude one another 
more strongly, and one is much less likely to find two electrons 
in a given small volume of jellium than in the same volume of 
spherium. 

We conclude from these comparisons that (at least) two 
parameters are required to characterize a uniform electron gas. 
Although the parameter choice is not unique, we believe that 
the first should be a one-electron quantity, such as the density 
p (or, equivalently, the Seitz radius r s ) and the second should 
be a two-electron quantity such as 77 = h(r,r), where h is the 
pair correlation function defined by 

p 2 (ri,r 2 ) = ip(r 1 )p(r 2 )[l+7j(n ! r 2 )] (28) 

and P2 is the diagonal part of the second-order spinless density 
matrix (TJ. 
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Accordingly, we propose to embark on a comprehensive 
study of L-spherium and L-glomium atoms, in order eventually 
to liberate the LDA from jellium's yoke through a process 
of radical generalization. The results of these spherium and 
glomium calculations will generalize the known properties 
of jellium, because we have shown that the energies of L- 
spherium and L-glomium approach those of 2D jellium and 
3D jellium, as L becomes large. 

In the remaining sections, we will confine our attention to 
the (two-dimensional) L-spherium atoms. However, exactly 
the same approach can and will be used to address the (three- 
dimensional) L-glomium atoms in the future. 



FIG. 2. Specific Coulomb holes for 0-spherium (dotted) with R = yl 
and 2D jellium (solid). Both are uniform gases with p = 1 /(14?r). 



V. LESSONS FROM SPHERIUM AND GLOMIUM 
A. A modest proposal 

The discovery that uniform electron gases are characterized 
by two parameters (p and 77) has many ramifications but one 
of the most obvious is that the foundations of the venerable 
Local Density Approximation need to be rebuilt. 

The traditional LDA writes the correlation energy of a molec- 
ular system as 



?ks , 



C(p)dr 



(29) 



thereby assuming that the contribution from each point r de- 
pends only on the one-electron density p (r) at that point. How- 
ever, now that we know that the energy of a uniform electron 



gas depends on p and 77, it is natural to replace Eq. (29 1 by the 
generalized expression 
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(30) 



In a sense, this two-parameter LDA represents a convergence in 
the evolution of Density Functional Theory (which stresses the 
one-electron density) and Intracule Functional Theory (which 
focuses on the two-electron density). 

How can we find the new density functional C(p,r\)l One 
could take an empirical approach but that is an overused option 
within the DFT community ll73l and we feel that it is more 
satisfactory to derive it from the uniform electron gas. As we 



show in Section V C it turns out that it is easy to compute 
the exact values of Ls, Ey, Ej and E\ for any L-spherium or 
L-glomium atom and, therefore, if one knew the exact wave- 
functions and energies of L-spherium and L-glomium for a 
wide range of L and R, one could extract the exact Kohn-Sham 
correlation energies 

eP=E-T s -E w -E s -E x (31) 

and determine the exact dependence of these on p and rj . 



B. Basis sets and integrals 

The Hamiltonian for L-spherium is 

in n 1 
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(32) 



and the natural basis functions for HF and correlated calcula- 
tions on this are the spherical harmonics Y/„,(0,0) introduced 
in Section|II] These functions are orthonormal ||64]| 

(Yim\Yi' m i) = 8/ ji 8 mjn t (33) 

and are eigenfunctions of the Laplacian, so that 

(Y lm \V 2 \Y,, m ,) = -1(1 + 1) 8 U , 8 mW (34) 

The required two-electron repulsion integrals can be found 
using the standard methods of two-electron integral theory 
J74). For example, the spherical harmonic resolution of the 
Coulomb operator I75if79ll 
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yields the general formula 
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where the one-electron integrals over three spherical harmonics 
involve Wigner 3 j symbols l64l and the sum over / and m is 
limited by the Clebsch-Gordan selection rules. 



C. Hartree-Fock calculations 

Unlike our calculations for electrons in a cube ll80l . the HF 
calculations for our present systems are trivial. Because the 
shells in L-spherium are filled, the restricted[? ] HF and KS 
orbitals are identical and are simply the spherical harmonics. 
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TABLE V. Reduced energies, densities and T) values of L-spherium with the four smallest eigenradii' R 



Wavefunction-based energies Kohn-Sham energies Ingredients of the new model 
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orbital-based correlation methods may be particularly effective. 
We now consider some of these. 

a. Configuration interaction H81]I This was the original 
scheme for proceeding beyond the HF approximation. It has 
fallen out of favor with many quantum chemists, because its 
size-inconsistency and size-inextensivity seriously hamper its 
efficacy for computing the energetics of chemical reactions. 
Nevertheless, for calculations of the energy of L-spherium, 
its variational character, systematic improvability and lack of 
convergence issues make it an attractive option. 

b. Mpller-Plesset perturbation theory [82 ] In an earlier 
paper we showed that the MP2, MP3, MP4 and MP5 
energies of 0-spherium can be found in closed form, for any 
value of R. However, although we observed that the MP series 
seems to converge rapidly for small R, its convergence was 
much less satisfactory for R > 1, where p < 0.15. Unfortu- 
nately, this is useless for our present purposes, because many 
of the L-spherium systems in Table |V| have much larger radii 
and lower densities than these values. 

c. Coupled-cluster theory R83\l The coupled-cluster hier- 
archy (viz. CCSD, CCSDT, CCSDTQ, . . . ) probably converges 
much better than the MPn series, and this should certainly be 
explored in the future, but we suspect that it will nonetheless 
perform poorly in the large-/?, small-p systems where static 
correlation dominates dynamic correlation |49 . 84, 85 1 and the 
single-configuration HF wavefunction is an inadequate starting 
point. 

d. Explicitly-correlated methods 4561/ The CI, MP and 
CC approaches expand the exact wavefunction as a linear com- 



These orbitals yield the reduced energy contributions 



Ts 



L(L + 2) 



Ej = +- 



Ek = 



L+l 
2R 



4L3 



4R 2 

L v = +0 

(L+l) 2 
R 

L, -1/2, 1/2, L + 2 
-L-l/2, 2, L + 3/2 
4(L+1) 



Lx = 



3kR 



(37) 
(38) 

(39) 
(40) 

(41) 



where 4L3 is a balanced hypergeometric function [64] of unit 
argument[? ]. It is encouraging to discover that Eq. (40 1 
approaches Eq. (41 1 in the large-L limit l57l . 

m 



We have used Eqs ( |37[ l - ( |40| l to compute the exact HF 
energies of L-spherium for a number of the eigenradii discussed 



in Section II C These, together with the Ts, Ej and Lx values, 



and p values from Eq. p), are shown in Table [V] 



D. Orbital-based correlation methods 



Although it is easy to find the HF energy of L-spherium, the 
calculation of its exact energy is not a trivial matter. How can 
this best be achieved? The fact that the occupied and virtual 
orbitals are simple functions (spherical harmonics), so the AO 
— >• MO integral transformation is unnecessary, suggests that 
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bination of determinants and it has been known since the early 
days of quantum mechanics that this ansatz struggles to de- 
scribe the interelectronic cusps [72|. The R12 methods over- 
come this deficiency by explicitly including terms that are 
linear in r, ; - in the wavefunction, and converge much more 
rapidly as the basis set is enlarged but, unfortunately, they also 
require a number of non-standard integrals. If these integrals 
are not computed exactly, they are approximated by resolu- 
tions in an auxiliary basis set. In either case, the computer 
implementation is complicated. 



E. Other correlation methods 

It is possible that the unusually high symmetry of L- 
spherium makes it well suited to correlation methods that are 
not based on the HF orbitals. We now consider two of these. 

e. Iterative Complement Interaction (ICI) method H87\ \88\l 
In this approach, the Schrodinger equation itself is used to 
generate a large set of «-electron functions that are then linearly 
combined to approximate the true wavefunction. It has been 
spectacularly successful in systems with a small number of 
electrons, but has not yet been applied to systems as large as 2- 
or 3-spherium (with 18 and 32 electrons, respectively). 

/ Quantum Monte Carlo methods [89. 90 J Of the var- 
ious methods in this family, Diffusion Monte Carlo (DMC) 
usually yields the greatest accuracy. In this approach, the 
Schrodinger equation is transformed into a diffusion equation 
in imaginary time T and, in the limit as T — > °°, the ground state 
energy is approached. Unfortunately, to be practically feasible, 
the method normally requires that the wavefunction's nodes be 
known and, despite some recent progress lr9"Tl , the node prob- 
lem remains unsolved. The quality of the trial wavefunction 
and finite-size errors are other potential restrictions ||92ll93l . 



F. Results and holes 

The discovery [51 1 that the Schrodinger equation for 0- 
spherium (and 0-glomium) is exactly solvable for each of its 
eigenradii is extremely helpful, for it allows us to generate the 
exact energies E and resulting Kohn-Sham correlation energies 
Eq s for the first four atoms in Table [v] without needing to 
perform any of the correlated calculations described above. 
However, these are the easiest cases, the "low-hanging fruit" 
so to speak, and there remain large gaps in the Table. We could 
have filled these gaps with rough estimates of the exact ener- 
gies but we prefer to leave them empty, to emphasize that there 
is much to do in this field and to challenge the correlation ex- 
perts in the wavefunction and density functional communities 
to address these beautifully simple systems. 

Once this is done, the data in the final three columns of Table 
|Vj will provide the ingredients for the construction of a new 
correlation density functional Eq S (p, tj) that will be exact for 
all L-spherium atoms and for 2D jellium. 

Finally, of course, an analogous strategy will yield a new 
functional that is exact for all L-glomium atoms and for 3D 
jellium. We will recommend that these functionals replace the 
LDA correlation functionals that are now being used. 

VI. CONCLUDING REMARK 

All uniform electron gases are equal, but some are more 
equal than others. 
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